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The colour fields, created by a static gluon-quark-antiquark system, are computed in quenched 
SU(3) lattice QCD, in a 24^ x 48 lattice at /3 = 6.2 and a = 0.07261(85) fm. We compute the hybrid 
Wilson Loop including the cases when the gluon and the antiquark are superposed, i. e., the quark- 
antiquark case and when the quark and antiquark are superposed, i. e., the gluon-gluon case. The 
Casimir scaling is investigated, in the two gluon glueball case the Casimir scaling is consistent with 
the formation of an adjoint string. Measuring the decay of the tail in the mid section of the flux tube 
for the two gluon glueball and for the quark- antiquark meson, we determined the penetration length 
and present a gauge invariant effective dual gluon mass of 0.905 =b 0.163 GeV. The coherence length 
was also studied using the dual Ginzburg-Landau approach. With the penetration and coherence 
lengths we determined the Ginzburg-Landau dimensionless parameter, this result is consistent with 
a type II superconductor picture. 



I. INTRODUCTION 

Understanding how the confinement arises from QCD 
is a central problem of strong interaction physics. In 
1970's, Nambu [Ij, 't Hooft [2j and Mandelstam [3j pro- 
posed that quark confinement would be physically in- 
terpreted using the dual version of the superconduc- 
tivity, the QCD vacuum state to behave like a mag- 
netic superconductor. In the ordinary superconductor, 
Cooper-pair condensation leads to the Meissner effect, 
and the magnetic flux is excluded or squeezed in a quasi- 
one-dimensional tube, the Abrikosov vortex, where the 
magnetic flux is quantized topologically. This confine- 
ment mechanism has indeed been established in compact 
QED [4-6J. From the Regge trajectory of hadrons and 
the lattice QCD results, the confinement force between 
the color-electric charge is characterized by the universal 
physical quantity of the string tension, and it is brought 
by one-dimensional squeezing of the color-electric flux in 
the QCD vacuum. Hence, the QCD vacuum can be re- 
garded as a dual version of a superconductor, based on 
the above similarities on the low-dimensionalization of 
the quantized flux between charges. 

In the dual-superconductor picture for the QCD vac- 
uum, the squeezing of the color-electric flux between 
quarks is realized by the dual Meissner effect, as the re- 
sult of condensation of color-magnetic monopoles, which 
is the dual version of the electrical charged Cooper pair. 

There is also evidence for the dual superconductor 
picture from numerical lattice QCD simulations [7]. If 
magnetic monopoles are condensed in the vacuum, then 
the electric sources are confined by electric flux tubes, 
as magnetic charges would be confined by Abrikosov- 
Nielsen-Olesen (ANO) vortices [8-lOj in an ordinary su- 
perconductor (Meissner effect). This is supported by ex- 
perimental observations like Regge trajectories, [llj, and 
lattice simulations showing a linear quark- antiquark po- 
tential. Therefore, the flux tubes correspond to vortices 
and the flux between charges must be confined. The 
confinement of magnetic field lines in superconductors 
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gluon mass in the literature. 



occurs through an Anderson-Higgs mechanism such that 
the photons of the electromagnetic field acquire a mass 
and are therefore exponentially damped in the supercon- 
ductor. In superconductors, the penetration length, A^, 
of the field in the London equation has a direct relation 
with an effective mass of the interaction particle fields, i. 
e., the photon. 

Thus, by measuring the penetration length of the fields 
in the QCD vacuum, we can estimate the mass of the dual 
gluon, if we further explore the analogy between QCD 
and superconductors. The dual gluon has been studied, 
[T2H19], however none of this studies have presented a 
value for the gauge invariant dual gluon mass in SU(3) 
lattice QCD. In Table [T| and [TT| we present some results 
found in the literature for the effective dual gluon mass 
and for the effective gluon mass, respectively. 

Bicudo, Cardoso and Oliveira, \45\ [46 j, established an 
analogy between the static potential and a type-II su- 
perconductor for the confinement in QCD, illustrated in 
Fig. [ij with repulsion of the fundamental strings and 
with the string tension of the first topological excitation 
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of the string (the adjoint string) larger than the double 
of the fundamental string tension. In type-I supercon- 
ductor the fundamental strings would be attracted and 
would fuse into an adjoint string. Cardoso et al., |47|[48] 
presented the color fields for this system and the results 
are compatible with the formation of an adjoint string 
between two gluon glueball. 

In this paper, we present the first value for the dual 
gluon mass in a SU(3) gauge independent lattice QCD 
and a detailed study of the Casimir scaling. In section 
II, we introduce the lattice QCD formulation. We briefly 
review the Wilson loop for this system, which was used 
in Bicudo et al. [46j, Cardoso et al. [45] and Cardoso 
et al. [48j, and show how we compute the color fields 
and as well as the lagrangian and energy densities distri- 
butions. In section III, the numerical results are shown. 
We present the results for the color field profiles in the 
mid flux tube section for the static hybrid gluon-quark- 
antiquark system, in a U shape geometry. A detailed 
study of the Casimir scaling is done. Measuring the 
penetration length of the tail in the mid flux tube, we 
present a value for the effective dual gluon mass. The co- 
herence length and the Ginzburg-Landau dimensionless 



parameter are studied using the dual Ginzburg-Landau 
approach. Finally, we present the conclusion in section 
IV. 
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Figure 1: String attraction and fusion, and string 
repulsion, respectively in type I and II superconductors. 



II. THE WILSON LOOPS AND COLOUR 
FIELDS 

The Wilson loop for the static hybrid gluon-quark- 
antiquark was deducted in [iF'iSJ, therefore we only 
present the fundamental expressions. The Wilson loop 
for this system is given by 



W1W2 - \W2^ 



(1) 



where Wi , W2 and are the simple Wilson loops shown 
in Fig. [2] Importantly, for the study of the Casimir Scal- 
ing, when n = 0, Wi = 3 and W2 = W3, the operator 
reduces to the mesonic Wilson loop and when ja = u and 
ri = r2 = r, W2 = Wl and W3 = 3, Wgqq reduces to 
Wgqq{r,r,t) = |W(r, t)p — 1, that is the Wilson loop in 
the adjoint representation used to compute the potential 
between two static gluons. 

In order to improve the signal to noise ratio of the 
Wilson loop, we use the APE smearing defined by 



+^^/7,(5)[/^(5 + z.)[/t(5 + /i)) , (2) 

with w = 0.2 and iterate this procedure 25 times in the 
spatial direction. 

To achieve better accuracy in the flux tube, we ap- 
ply the hypercubic blocking (HYP) in the time direction, 
[49], with 



ai = 0.75, a2 = 0.6, 0^3 = 0.3 



(3) 



We obtain the chromoelectric and chromomagnetic 
fields on the lattice, by using. 



(Ef) = (Poi) 



{W) 



(4) 
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Figure 2: (a) Wilson loop for the gqq and equival ent 
position of the static antiquark, gluon, and quark, (b) 
Simple Wilson loops that make the gqq Wilson loop. 



and, 



/b2X _ jWPjk) ,p > 



(5) 



where the jk indices of the plaquette complement the 
index i of the magnetic field, and where the plaquette is 
given by 

P^, (s) = 1 - ^ReTr [U^{s)UAs + ^i)Ul{s + v)Ul{s)\ . 

(6) 

The energy (H) and lagrangian (C) densities are given 

by 



(7) 



(8) 



Notice that we only apply the smearing technique to the 
Wilson loop. 



III. RESULTS 

Here we present the results of our simulations with 
286 24^ X 48, /3 = 6.2 quenched configurations generated 
with the version 6 of the MILC code [50j , via a combina- 
tion of Cabbibo-Mariani and overrelaxed updates. The 
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(a) U shape geometry. 



(b) L shape geometry. 



Figure 3: gluon- quark- antiquark geometries, U and L 
shapes. 



results are presented in lattice spacing units of a, with 
a = 0.07261(85) fm or a'^ = 2718 ± 32 MeV. 

In this work two geometries for the hybrid system, 
gluon-quark-antiquark, are investigated: a U shape and 
a L shape geometry, both defined in Fig. [3j Note that, in 
the L shape geometry only one case is studied here, the 
case when the gluon and the antiquark are superposed, 
ri = 0, this case corresponds to the quark- antiquark case. 
This particular case is relevant to study the Casimir scal- 
ing. 

The use of the APE (in space) and HYP (in time) 
smearing allows us to have better results for the flux tube 
and reduce large plaquette fluctuations. Notice in Fig. [4] 
the improvement of the signal in the tail of the flux tube 
profile. 



A. U profiles and Casimir scaling 

In Fig. [5] we present the profiles for the U geometry 
with / = 8 and d between and 12 at ^ = 4. We can 
see the stretching and partial splitting of the flux tube 
in the equatorial plane {y = 4) between the quark and 
antiquark and z = 0. When (i = 0, the quark and the 
antiquark are superposed and this corresponds to the two 
gluon glueball case. 

For d = 2 and 4 at = 4, the separation between the 
two flux tube is not visible, this is due to the overlap in 
the tails of the flux tube which contributes for the total 
field and, for large separations, the tails of the flux tubes 
contribute to a non zero field at x = 0. 

We measure the quotient between the energy density of 
the two gluon glueball system and of the meson system, 
in the mediatrix plane between the two particles {x = 
0). The results are shown in Fig. |6] In Fig. [6a| we 
present the results for r = at x = z = and in Fig. 
[6b| we present the results for r = (x, z) at = 4. We 
make a constant fit to the data in Fig. |6j the result 
for Fig. [6a| is 2.25096 ± 0.0244972 and for Fig. |6b] is 
2.23591 ± 0.0598732. As it can be seen, these results 
are consistent with the Casimir scaling, with a factor 
of 9/4 between the energy density in the glueball and 
in the meson. This corresponds to the formation of an 
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Figure 4: Comparison between results for lagrangian 

density with APE smearing only in space and for 
lagrangian density with APE smearing in space and 

Hypercubic Blocking (HYP) in Time for the U 
geometry with d = and / = 8. In |(a)[ the lines were 
drawn for convenience and therefore do not represent 
results from any kind of interpolation. 



adjoint string between the two gluons, and in the quark- 
antiquark case we have a fundamental string. The points 
outside the field, i.e., as the field approaches to zero, were 
discarded. 
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Table III: Results for the dual gluon mass, where (la) is 
for the two gluon glueball and (2a) for the 
quark- ant iqu ark cases at y = 4 and z = with 
r = X > 2 and z = 0, and (lb) is for the two gluon 
glueball and (2b) for the quark-antiquark cases ai y = 4 
with r = (x, z) > 2. 



tial rather than the gaussian. This result reinforces the 
correspondence between the QCD vacuum and the Dual 
Superconductor. Baker et al. [52] presented as solution 
to the electric and magnetic colour flux tube the modi- 
fied Bessel functions. We also fitted our data with the 
modified Bessel function of zero order, Kq {iir). 

So in our case, the mean value of the squared fields 
should decay as e~^^^ or like Kq {jir)^ where /i should 
be the effective dual gluon mass. However, since we have 
a flux tube and not a wall, i. e., the case study is not 
planar, the correct fit is the modified Bessel function of 
zero order, suitable for cylindrical flux tubes. 

Since ^ and C are all gauge invariant, 

we can compute a value for the effective dual gluon mass 
that is gauge independent by studying the decay of one of 
these quantities in the middle of the flux tube, between 
the sources. Fitting the chromoelectric field and the la- 
grangian density section in the mid distance of the flux 
tube of the meson and the two gluon glueball, we obtain 



the results presented in Table III for the effective dual 
gluon mass. The results are compatible with an identical 
decay of the Lagrangian density in both cases and are of 
the order of ~ 1 GeV. 



In common superconductivity the magnetic field de- 
cays with B ~ ^-"^1^1^ and this could be interpreted in 
terms of an effective mass for the photon = l/A^. 
Some studies have pointed a similar behavior in QCD, 
[13, 51-53J. as in a superconductor, the chromoelectric 
field decays as ~ ^-"^l^i^ as they depart from the singular 
vortex flux tube, where r stands for the distance away 
from the center of the flux tube. This is consistent with 
the dual superconductor picture. 

We measure the decay of the mid flux tube section in 
the quark-antiquark and the two gluon glueball. We tried 
to fit our data by a gaussian and by a exponential, but 
the /dof clearly indicated that the decay is exponen- 



C. Ginzburg Landau Equations 

The failure of fitting the entire set of the chromoelectric 
in the mid flux tube data hints at the existence of another 
scale, in addition to the dual gluon mass, /i = A~^. The 
extraction of this second scale, the coherence length, ^, 
can be done using the effective dual Ginzburg-Landau 
approach. The ratio between the penetration length and 
the coherence length. 
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Figure 5: Results for the U geometry at y = 4 and z = 0. 
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yure 6: Results for the two gluon glueball {d = and / = 8, U geometry) energy density over the meson (ri 
and r2 = 8, L geometry) energy density at x = and Casimir scaling factor of 9/4 (broken line). 
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Table IV: Results obtained by fitting the numeric solution of the Ginzburg-Landau and Ampere equations to the 
colour fields and lagrangian density data with A fixed, gg - two gluon glueball results, qq - quark- ant iquark results in 

XZ plane at y = 4. 



is the famous Ginzburg-Landau dimensionless parameter 
used to determine whether a superconductor is type I or 
type 11. 

The Ginzburg Landau equation and the Ampere equa- 
tion are given by 

^(V - ^A)2^ ± - = , (10) 

2m he 



Vx Vx A= :^(^*V^- W^*-^|^pA). (11) 



2imc 



he 



Defining a = ^ c2 ^ A 



a ' ^ 



^ and a 



-^A, we arrive at the simplified equations 



^^(V ± za) V + - ^IV^I V = 



(12) 



A^'V X V X a = 



2i 



iV^pa . (13) 



In the case of a vortex, we have, in cylinder coordinates 
a(r) = a^{p)e^ and ?/^(r) = ii(/))e~*'^, so the equations 
are given by 

e{-^{pp)--,+2^u-alu)=u'-u, (14) 
^ pap dp p ^ ' 



Id du 
pdp^^dp' 



1 2.1 
^-p 



(15) 



Solving these equations with A fixed, obtained from 
the modified Bessel function of zero order fit in Table 
Jnj and fitting the numerical result to the the results of 
the colour fields and the lagrangian density in the middle 
of the fiux tube, in XZ plane, we obtain the results in 
Table |IV| and Fig. [7[ Therefore, the Ginzburg-Landau 
parameter, is greater than 2~^/^, which corresponds 
to a type II superconductor. 



IV. CONCLUSIONS 

We present the first gauge invariant result for the effec- 
tive dual gluon mass in a pure gauge SU(3) QCD lattice. 
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(a) Two gluon glueball case. 
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(b) Quark-antiquark case. 

Figure 7: Fit to the chromoelectric field data using the 
numeric solution of the Ginzburg-Landau and Ampere 
equations, red points. The blue points are the results 
for chromoelectric field. 



A detailed study of the Casimir scaling is presented and 
the results are consistent to the formation of an adjoint 
string between the two gluon and agrees with the Casimir 
Scaling measured by Bali [54J, with a factor of 9/4. 

We also study the decay of the fields perpendicular to 
the flux tube, and obtain an exponential decay consis- 
tent with the dual superconductor model of confinement. 
By considering the decay factor of the mid distance tube 



7 



flux, we obtain a gauge invariant result of ~ 1 GeV for 
the effective dual gluon mass and this value is consistent 
with previous analytical theoretical calculations and in- 
dependent phenomenological estimates for the effective 
dual gluon mass and gluon mass found in the literature. 

Applying the Ginzburg Landau and the Ampere equa- 
tions, we are able to perform a complete fit to the colour 
fields in the mid flux tube and determine the coherence 
length parameter, with the penetration length fixed and 
determined from the modified Bessel function of zero or- 
der fit. Finally, with this two parameters, we determine 
the Ginzburg-Landau dimensionless parameter, which 



indicates a type II superconductor region see Table IV 



Nevertheless, SU(3) QCD is not completely equivalent to 
a dual superconductor. 

We find that all six components Ef and Bf have similar 
flux tube profiles, while in a superconductor only the 
component of the magnetic field B parallel to the flux 
tube has such a profile. The dominant component of 
the chromoelectric field, E^^ parallel to the flux tube. 



is about two times the E^ and E'^ components. In the 
chromomagnetic field the B'^ and B^ are dominant and 
are about is about 1.5 times the By component in the flux 
tube. The dominant component of the chromoelectric 
field is about 3.3 times the ^nd 2.2 the \Bl \ and 
in this region. This suggests in QCD the existence of an 
gluon mass, identical to the dual gluon mass. Combining 
all Ef and Bf components we find the result for the dual 
gluon mass and gluon mass of 0.905 ± 0.163 GeV. 
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